Fluctuator model of memory dip in hopping insulators 
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We show that the non-equilibrium dynamic in two-dimensional electron glasses close to metal- 
dielectric transition is dramatically sensitive to electric fields confinement inside the sample, which 
leads to a nearly thermally activated conductance behavior and a strong non-equilibrium conduc- 
tance response to the gate voltage, i. e. memory dip in a field dependence of conductance. Theory 
accounts qualitatively and quantitatively for the universal temperature and field dependencies of 
memory dip. 

PACS numbers: 73.23.-b 72.70.+m 71.55.Jv 73.61. Jc 73.50.-h 73.50.Td 



Slow relaxation in glassy materials is associated with 
transitions between local minima in configurational space ® 

separated by high potential barriers. At low temperature ea^< ® 

a glassy material does not approach the global energy 
minimum, but reaches some local minimum. There exist 
other local minima possessing lower energies, where the 
system transfers after certain expectation time. Reason- 
ably broad distribution of transition barrier heights and 
lengths results in the logarithmic character of relaxation 
because transition rates depend on these parameters ex- 
ponentially. 

At low temperature (4 — 20K and below) this loga- 
rithmic relaxation is observed in a variety of disordered 
materials. Particularly, a sudden application of an ex- 
ternal electric field to amorphous dielectrics results in a 
fast increase of dielectric constant with its subsequent 
logarithmic relaxation back to equilibrium [l|, |2|. Sim- 
ilar behaviour of conductance is observed after sudden 
application of a gate voltage in two-dimensional Ander- 
son insulators including indium oxide [3-5], granular alu- 
minium [a-lSj and ultrathin films of Bi and Pb. [9| 

A general mechanism for the non-equilibrium relax- 
ation of conductance, earlier used to interpret experi- 
mental data for dielectric constant in glasses [3, llOJ , can 



FIG. 1: Reduction of DOS due to the interaction of electron 
(blue circle with — ) and configurational transition (fiuctuator, 
shown by the black arrow with a characteristic dipole moment 
ea) interaction. If the interaction of fluctuator and dipole is 
very large escape of electron (blue arrow), must be accompa- 
nied by the simultaneous transition of fluctuator (state I to 
state II). 



trie constant. Thus equilibration of system after its dis- 
turbance, e. g. by the gate voltage, results in slow reduc- 
tion of conductance by means of trapping of low energy 
electrons, primarily released by the gate voltage. 

Another interesting observation related to glassy prop- 
erties is the memory dip seen in conductance after a fast 
sweep of the gate voltage around its equilibrium value, 
V^ = (see l^-shaped conductance dependence on that 
voltage in Fig. [5]). The greater the voltage the stronger 
the system departs from equilibrium increasing the con- 
ductance, according to the previously described scenario. 

(1/2) 



be described as following [ITl, [13 ■ Application of external Memory dip can be characterized by its half- width V^ 



field (gate voltage) reduces system stability in the present 
local minimum and creates new pathways for its configu- 
rational transitions to lower energy states. Reduction of 
system energy due to configurational transitions lowers 
the electronic density of states (DOS) and consequently 
system conductance. 

DOS decreases due to an interaction of configurational 
transitions with electrons in localized states (see interac- 
tion of electron with local transition shown in Fig. [ij. 
If this interaction is sufficiently strong then the electron 
needs the large energy of that interaction to leave or enter 
this state. At low temperature such process becomes for- 
bidden and the only possible excitation of electron must 
be accompanied by the configurational transition. Since 
configurational transitions occurs very rarely this elec- 
tron cannot efficiently respond to the external electric 
field, i. e., contribute to hopping conductivity or dielec- 



and depth 5Cmax/C . 

The temperature dependence of conductance, C, mem- 
ory dip parameters, Vg \ and SCmax/G, and logarith- 
mic relaxation rate of conductance, r/„ = dln(C)/dln(i), 
contain significant information about the nature of the 
glassy state. It turns out that within the temperature 
range of interest, 4K< T < 25K, the conductance of 
most of samples shows nearly thermally activated behav- 
ior, C ex e-(^o/T)"^ 0.8 < n < 1 0,1,11]. This behaviour 
can be interpreted as a consequence of field confinement 
in 2 dimensions {2D) 13, |lj|, rather than the standard 
variable range hopping mechanis m, 1/4 < n < 1/2, con- 
sidered in earlier work llj, |l2|, Il6l421|. In this letter 
we suggest theoretical model of electronic glass involving 
the 2D field confinement, which is capable to interpret 
several observed behaviors at intermediate temperature, 
T > 4:K, including (A) linear temperature dependence 



of a memory dip half- width and its universahty with re- 
spect to a sample conductance [8| (cf. [5|); (B) nearly 
T^^ dependence of AC/C varying for samples with dif- 
ferent conductances Q; and (C) nearly T~'^l'^ tempera- 
ture dependence of logarithmic relaxation rate [9|]- Ab- 
solute values of these parameters agree with the common 
sense expectations. 

If the dielectric constant of the film of interest, Ki„ ex- 
ceeds that of the environment, kj„ 3> n^xx = ('«i +K2)/2 
(ki, K2 are dielectric constants of surrounding materi- 
als), the interaction of electrons differs from the standard 
Coulomb interaction, e^ /{nir,r), because of the field con- 
finement within the layer [IJ, [ij] (see also Ref. [15| . 
where the field confinement was considered for Ti oxide 
based materials) . In that case the interaction of two elec- 
trons at short distances r < d/2 can be expressed with 
the logarithmic accuracy as 



Uc(r) 
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This interaction results in the hard Coulomb gap A* in 



a DOS leading to the conductance behaviour |13lll4| 



a = aoe 



En = 2.i 



(2) 



where a stands for the electron localization length. 

This dependence is similar to the one reported in Ref. 
[8| for various granular aluminium samples in the tem- 
perature domain of interest 4K < T < 20K and the most 
significant temperature dependence in exponent is associ- 
ated with the Arrhenius term in the exponent A^/{kBT). 
Consequently, one can expect that the non-equilibrium 
conductance behavior is also associated with this term. 
Then the relaxation and memory dip are determined by 
the dielectric constant time and field dependencies as 



da{Vg,t) _ A, SK,niVg,t) 

a ksT Kin 



(3) 



It is quite natural to expect that the dielectric constant 
of Anderson insulators behaves similarly to that in amor- 
phous solids [l|, y, llO| yielding experimentally observed 
conductance behavior. 

To characterize the non-equilibrium dynamics one has 
to describe configurational transitions. In dielectric 
glasses relaxation and ageing phenomena are associated 
with the two-level systems (TLS) [24I formed by tunneling 
transitions of atoms or groups of atoms between two close 
energy minima separated by potential barriers. TLS pos- 
sess a universal statistics with respect to their energies, 
E, and relaxation times, r, 



P{E,r) = ^. 
and a small dipole moments /istr ~ 1 — 3D. 



(4) 




Vg(V) 



FIG. 2: Typical shape of memory dip 



One can expect that in Anderson insulators electrons 
themselves can form the electron glass state 18|, l20|, l21 , 
l2S - t25j |. Assuming that electronic glass is formed due to 
Coulomb interaction one can estimate the glass transi- 
tion temperature, Tg is given by ksTg = Eq (see Eq. ([2]), 
cf. Ref. [20|). In this state the slow dynamics is associ- 
ated with collective electronic transitions between close 
energy minima. The logarithmic relaxation can be due to 
exponential sensitivity of electron transition rate to the 
hopping length and the number of participating electrons 
[12I llGJ I . Transitions observable within the time domain 
of the typical experiment, i. e., from several seconds to 
several weeks, cannot involve more than 5 — 10 electrons 
|16| located close to each other. Then we can still use a 
TLS like model Eq. (j4]) to characterize electronic tran- 
sitions as local fluctuators [16| possessing density, P and 
dipole moments /i* ~ ea; the localization length a is the 
only length parameter of the problem. It is convenient 
to characterize fuctuators by the dimensionless product 
P*fJ-l/Kin « x*i so one has 



ea, P, 



X*^*' 



(5) 



In amorphous solids x* is a universal parameter of order 
of 10^"^ — 10"'', which can possess weak (logarithmic) 
temperature dependence. As we will see below, the sim- 
ilar assumption leads to the right estimate of the mem- 
ory dip depth in conductance. Although the mechanism 
of electron glass formation is beyond the scope of the 
present work we believe that the long range 3 — D inter- 
action is important there (cf. [2|). Therefore the forma- 
tion of glassy state can be much less efficient in the case 
of small concentration of carriers (less than one per cu- 
bic thickness of material) in agreement with experimen- 
tal observations [5|. Also if localization length exceeds 
the thickness the effective interaction is two-dimensional, 
which can also suppress the glass formation [16|. 

Consider the effect of fiuctuators on the dielectric con- 
stant determining the non-equilibrium behavior Eq. ([3|). 
The dielectric constant is associated with the polarization 
of electron-hole pairs of various lengths r, smaller than 
the hopping length. Similarly to Ref. [10[ the part of di- 



electric constant most sensitive to its interaction with 
fluctuators is associated with low-energy electron-hole 
pairs separated by intermediate distance, r, such that 
a < r < d/2. It is harder to disturb shorter dipoles due 
to their large tunneling splitting, Aq > fc^To » fcsT, 
while longer pairs do not make important contribution 
due to the "hard" Coulomb gap at corresponding ener- 
gies caused by field confinement in 2D. Let those dipoles 
be characterized by the density function F(A, r) depend- 
ing on their energy, A and size r. This size determines 
the tunneling amplitude of electron between two local- 
ized states, Aq « EQe~^^°'. Then one can approximate 
the contribution of interest to the dielectric constant as 

(c. f. [ii,!!!) 
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Two contributions in the right hand side of Eq. ([6|) are 
identical to resonant and relaxational contributions to 
the dielectric constant in amorphous solids [2, llOj and 
the latter contribution is more significant because it is 
determined by lowest energies of pairs. 

In the absence of fluctuators the density function of 
sufficiently large dipoles, a < r < d\s given by [26|, [27 1 



F(A,r) 



107r2 



(7) 



The interaction of dipoles with fluctuators reduces their 
density of states [lO|, [l^ • Only interaction exceeding the 
thermal energy is important. We assume, that the dis- 
tance between the fluctuator and closest charge of the 
dipole is smaller than the size of the dipole r, which is 
necessary to make the interaction larger than the thermal 
energy. This is true for r ^ dl2 and T > 4K [3]. Then 
the correction to the dipole density can be expressed as 
nl [l3| (remember, E = ^ A^ + Al) 



(SF(A,r) 
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The time dependence is due to a logarithmically uniform 
distribution of fluctuator relaxation time (Eq. Q), t 
stands for the observation time and r,„ is some charac- 
teristic minimum time, determined by the rate of a gate 
voltage application. 

Substituting Eq. ([8|) into Eq. ([6]), performing the 
straightforward integrations and substituting the final re- 
sult into Eq. ([3]) we obtain the following expression for 
the fluctuator correction to the conductance 



5C _ 32V^x* A, 
C ^ 75 



fcsTV ksT \2a) \t 



(9) 



This equation describes the ageing effect as well as the 
conductance time dependence the application of a large 
gate voltage bringing all relevant fluctuators out of equi- 
librium. The temperature dependence of the conduc- 
tance logarithmic relaxation rate r — (iln(C)/(iln(i) is 
close to the power law r oc T^-^/^ that is consistent with 
the observations of Ref. [9!| in ultrathin films of Bi and 
Pb. 

Consider the shape of the memory dip (Fig. [T]). The 
energy change of fluctuator with the dipole moment /i as- 
sociated with the application of a gate voltage Vg is given 
by 5E — Fg^cos{6) sa VgKifiCos{9)/{Kindtot), where 6 is 
the angle between directions of the dipole moment and 
external electric field Fg, dtot is the distance between 
sample and gate electrode and ki is the dielectric con- 
stant of the insulating layer separating the gate electrode 
from the sample (see Refs. [5|,|8| for details). Only fluctu- 
ators with sufficiently small energy E, < E < SE, will 
be really disturbed from equilibrium. This condition sets 
the upper constraint SE for integration over energy in Eq. 
dS]). Also the time dependent logarithm \n{tmax/tmin) in 
Eq. ^ should be modified. It is determined by the con- 
tribution of fiuctuators removed from equilibrium by the 
gate voltage, which possess the relaxation time smaller 
than the maximum time imaai determined by the time 
the sample was kept at fixed temperature after cooling 
and larger than the minimum time, tmim determined by 
the gate voltage sweep rate, r^, as tmin ^ kBT/{ers). 
Then the memory dip can be described by the equation 
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(10) 



Using this result one can estimate the half width of the 
memory dip as 



T/(l/2) 



16 



e Kia 



(11) 



This result agrees with the experimentally observed 
linear temperature dependence of the memory dip width 
at temperatures T > 4K [a, |8| . At lower temperatures 
the dependence is getting stronger. In our opinion this 
is because the estimate for gate voltage induced electric 
field, Fg « Vg Ki/{Kindtot), is no longer valid. One 
can use that macroscopic expression if the density of elec- 
trons injected due to the application of a gate voltage ex- 
ceeds one electron per the cubic sample thickness. This 
is satisfied for granular Aluminium samples down to 4K 
6], but fails at lower temperatures, where this field is 
defined by the injected electron, closest to the given fluc- 
tuator. One can show that at low temperatures a width 
of memory dip is proportional to squared temperature. 



Vc 



(1/2) (fc^T)^ d'^diKl 



The quantitative universality of the memory dip 
halfwidth discovered in [8J for samples with conductances 
different by orders of magnitude requires understanding. 
It is important that in the vicinity of a metal-insulator 
transition the ratio of two diverging parameters, tin/ a, is 
expected to be approximately constant on the dielectric 
side both in accordance with the theoretical analysis and 
experimental data [28|. One can naively expect that at 
the scale of localization length the characteristic kinet- 
ics and Coulomb energies should be close to each other, 
Eq = e^ / {Kind) ~ It? /{raa?) which, results in the univer- 
sal ratio a/€in ^ h^ /{rae^) k, O.SA. This is simply the 
Bohr's radius, as {m is the effective mass of electron, 
which is taken to be equal to the bare electronic mass). 

Under these assumptions Eq. dTl]) results in the uni- 
versal dependence of the memory dip halfwidth on the 
gate voltage Vg = vfT with the proportionality coef- 
ficient coefficient 77 depending only on the thickness dtot 
and dielectric constant ki of the layer separating the sam- 
ple and the gate electrode. For granular Aluminium sam- 
ple one has ki « 9, dtot — lOOOA which yields rj — 0.29V 
K^^ in excellent agreement with the experimental result 
ry« 0.25V K-i Q. 

The same approach does not work so good for the half- 
width of the memory dip reported in Ref. [5| for In203_a: 
sample having a relatively small resistance of 200kr2. We 
believe that this is because the localization length in this 
almost conducting sample exceeds the thickness of the 
sample dino ~ 30A and therefore one should replace 
the localization length a in Eq. pT|) with the sample 
thickness. Assuming that the localization length exceeds 
the sample thickness by an order of magnitude one can 
make the experimental data consistent with the theory. 
Other data reported in Ref. J5| correspond to very low 
temperatures. However, in granular Aluminium one can 
estimate (fitting activation energy 20 — 40K with Eq. ([T]) 
and using the relationship a sa 0.5Ki„(A)) that IhA < 
a < 50 A in all samples of thickness 200A so our estimate 
is justified. 

The temperature dependence of the depth of the mem- 



ory dip in Eq. (|10|) in the limit of Vg — > 00 is determined 
by power law temperature dependence T~^'^, possible 
weak temperature dependence of the dimensionless pa- 
rameter Y, and the logarithmic factor In ( ^^-^^^J^'^"^ ) . 
Integration of these three dependences could be respon- 
sible for the T~^ behavior of the memory dip depth re- 
ported in Ref. [7[. 

According to the experiment [7| the absolute value of 
the relative correction reaches few percents at T ~ AK 
for less conducting sample, which is consistent with Eq. 
(ITUl) if one set x* '^-^ 5 • 10"'* in agreement with the typi- 
cal values of that parameter typical value in amorphous 
solids [2|. 

Thus we suggested the model of electronic glass. In 
this model slow dynamics can be characterized by fiuc- 



tuators similar to TLS in amorphous solids. Our model 
is consistent with existing experimental data both quali- 
tatively and quantitatively. Theoretical predictions, Eqs. 
([TU|. f dTTj) can be tested varying system parameters in- 
cluding electron localization length, sample thickness and 
temperature. 
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